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AN AXISYMMETRIC PISTON-CYLINDER ARRANGEMENT 


by 
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Lewis Research Center 
Cleveland, Ohio 

and 

Gene E. Smith and George S. Springer 
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Ann Arbor, Michigan 


ABSTRACT 

The implicit factored method of Beam and Warming was employed to describe 
the flow and the fuel-air mixing in an axisymmetric piston-cylinder configu- 
ration during the intake and compression strokes. The governing equations 
were established on the basis of laminar flow. The increased mixing due to 
turbulence was simulated by appropriately chosen effective transport proper- 
ties. Calculations were performed for single-component gases and for two- 
component gas mixtures. The flow field was calculated as functions of time 
and position for different geometries, piston speeds, intake-charge-to- 
residual-gas-pressure ratios, and specie mass fractions of the intake charge. 
Results are presented in graphical form which show the formation, growth, 
and break-up of those vortices which form during the intake stroke and the 
mixing of fuel and air throughout the intake and compression strokes. It is 
shown that at bore-to-stroke ratio of less than unity, the vortices may 
break-up during the intake stroke. It is also shown that vortices which do 
not break-up during the intake stroke coalesce during the compression stroke. 
The results generated were compared to existing numerical solutions and to 
available experimental data. 


I. INTRODUCTION 

An understanding of the flow patterns inside cylinders of spark ignited 
reciprocating engines is required in developing efficient, low emission 
engines. For this reason, numerous investigators have studied flows inside 
the cylinders of reciprocating engines by both experimental (e.g., see survey 
in ref. 1) and numerical (e.g., see surveys in refs. 1 to 4) techniques. 
However, owing to the complexity of the phenomena, a complete description of 
the flow cannot be obtained by presently available experimental and numerical 
methods. Instead, progress towards the solution must be made by focusing 
attention on certain aspects of the problem. In this study, the flow pat- 
terns and mixture distributions were investigated inside an axisymmetric 
piston-cylinder configuration during the intake and compression strokes. 

The results were generated by a numerical method of solution of the governing 
equations. 
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II. DESCRIPTION OF THE PROBLEM 


The following problem was analyzed. A hollow circular cylinder {inner 
radius rp) is closed on one end by a flat piston and on the other end by a 
flat plate (fig. 1 ). The piston is connected to a crank shaft (radius r^) 
through a connecting rod (length !(-)• The piston is driven by rotation of 
the crank shaft about the crank pin at an angular velocity of n, resulting 
in a piston velocity of Up. 

The flat plate has a centrally located annular opening (valve opening, 
inner radius r^, outer radius r^) in it which opens instantaneously at the 
beginning of the intake stroke (crank angle ^ = 0) and closes instantane- 
ously at the end of the intake stroke [i = it). 

The temperatures at the cylinder wall T^., valve T^, cylinder head T^, 
and piston Tp are constants, but may have different values. 

The fluid^enters the piston-cylinder configuration described above 
(henceforth referred to as the cylinder) through the valve opening during 
the intake stroke. The static temperature T-j and static pressure P-j of 
the entering fluid are both constant. At the valve opening (seat angle a), 
the entering fluid may have velocity components in the radial (Vj.) and axial 
(V7) directions, but not in the tangential direction. The magnitude of the 
velocity at the valve opening depends upon the instantaneous flow field in- 
side the cylinder as well as the static temperature and pressure of the fluid 
at the valve opening. 

The viscous and thermally conducting fluid which enters the cylinder (in- 
take charge) is either a single-component ideal gas or a two-component non- 
reacting ideal gas mixture. The two-components of the intake charge are to 
simulate air and fuel. The thermodynamic and transport properties of the 
intake charge (except for the binary diffusion coefficient) are taken to be 
constants. The binary diffusion coefficient between the two-components, D/\g, 
is taken to be inversely proportional to the local mixture density p. The 
mass fraction of the two-component ideal gas mixture (Xa and Xg) may vary 
arbitrarily with time at the valve opening. 

At the beginning of the intake stroke, the gas in the clearance volume 
(residual gas) is taken to be a stagnant, single-component ideal gas at sta- 
tic temperature T-j and static pressure Pq (Pq ± Pi). The residual gas has 
the same thermodynamic and transport properties as component A of the intake 
charge. 

The problem just described is symmetric with respect to the axis of the 
cylinder. Accordingly, the tangential velocity is zero everywhere in the 
cylinder, and the radial velocity and the first-order radial derivatives of 
density, mass fractions of components A and B, axial velocity, and energy are 
zero at the center line. 


III. METHOD OF SOLUTION 

The basic equations governing the problem are the conservation equations 
of mass, radial momentum, axial momentum, and energy and the specie balance 
equations (refs. 5 and 6). To render these equations more amenable to numer- 
ical methods of solution, the following simplifications were made: (a) pres- 

sure gradient contribution to specie diffusion and momentum transport due to 
specie diffusion were not considered; (b) irreversible coupling effects be- 
tween the temperature gradient and specie concentration gradient (Soret and 
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Dufour effects) were neglected; (c) body forces were neglected; (d) the bulk 
viscosity was taken to be zero; and (e) radiation heat transfer was neglected. 

The governing equations, applicable to laminar flows, are summarized in 
table 1. The increased mixing due to turbulence was simulated by appropri- 
ately chosen effective transport properties (Section IV). Equations (1) to 
(19) constitute a closed system in the five basic dependent variables: den- 

sity (p), mass fraction of component A (X/^), radial velocity (V^), axial vel- 
ocity (V^), and energy (e). 

The boundary and initial conditions corresponding to equations (1) to 
(19) are given in table 2. The no-slip condition requires the fluid velocity 
next to a solid wall to be equal to the velocity of the solid wall (eqs. (20) 
and (21)). 

The gas temperature next to solid walls equals the temperature of the 
walls. With all walls maintained at constant temperatures, this boundary 
condition results in equations (22) to (25). 

The cylinder walls are assumed to be nonporous. This boundary condition 
is expressed by equations (26) and (27). The static pressure and static tem- 
perature of the gas at the valve opening are constant with respect to time 
as reflected by equation (28). The mass fraction of component A at the valve 
opening is expressed by equation (29). The radial and axial velocities at 
the valve opening are related, as indicated by equation (31). The axial vel- 
ocity at the valve opening is determined by applying the conservation of mass 
equation at the valve opening (eq. (32)). 

Furthermore, the symmetry conditions at the center line result in equa- 
tion (33). At time t equal to zero, the piston is at the TDC position 
{i = 0) and the residual gas in the clearance volume is a stagnant single- 
component ideal gas at static pressure and static temperature T-j. 

These initial conditions are specified by equation (34). 

Solutions to the governing equations and the corresponding initial and 
boundary conditions formulated above must be obtained by numerical methods. 

In this investigation, the implicit-f actor ed, finite-difference method of 
Beam and Warming (refs. 7 and 8) was used to obtain solutions. The detailed 
steps of the method was described elsewhere (refs. 9 and 10) and will not be 
repeated here. Essentially, the continuous domain inside the cylinder was 
represented by a grid system. The grid size in the axial direction was uni- 
form but changed with time as the piston moved towards or away from the cyl- 
inder head. The grid spacing in the radial direction did not change with 
time, but varied in size to resolve regions where the radial derivatives were 
steep. The governing equations were first transformed to a moving coordinate 
system corresponding to this grid and were then expressed in finite-differ- 
ence form. The resulting system of finite-difference equations had block 
tridiagonal coefficient matrices and hence could be solved by the Thomas 
algorithm (refs. 11 and 12). In addition, the boundary conditions were ex- 
pressed in a form consistent with the implicit factored method of Beam and 
Warming and the Thomas algorithm. A complete description of the equations 
may be found in reference 9. A copy of the computer code may be obtained 
from the authors. 


IV. RESULTS 

Numerical solutions were obtained to explore the flow patterns and the 
mixture distributions inside the cylinder during the intake and compression 
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strokes. Two types of problems were studied: (a) single-component gas 

flows, and (b) two-component gas mixture flows. 

The ranges of the parameters for which numerical solutions were obtained 
are summarized in table 3. The values of the parameters describing the geom- 
etry and operating conditions were selected so as to correspond to those 
typically found in spark-ignition engines. In selecting thermodynamic and 
transport properties, two constraints were taken into account. First, the 
values were selected so as to be physically reasonable, as will be discussed 
subsequently. Second, the values were chosen so as to permit the use of con- 
venient grid sizes (ref. 9). 

For problems involving a single-component gas, the thermodynamic proper- 
ties of the gas were those of air at temperature T-j. For problems involving 
a two-component gas mixture, the thermodynamic properties were those of air 
at temperature T-j for component A, and those of octane at temperature T-j 
for component B. In every problem, T-j was the temperature of the intake 
charge and the initial temperature of the residual gas. 

The value of the viscosity, y, was based on the observation that the ef- 
fective viscosity for turbulent flows inside spark ignition engine cylinders 
is roughly 100 times higher than the viscosity corresponding to laminar flows 
(ref. 13). The laminar viscosity of air at temperature T-j (340 K) and atmo- 
spheric pressure is 3.203xl0“5 kg/m-s. The value of the effective viscosity 
used is 100 times this value (table 3). 

The values of the mass diffusivity, D^Bj and thermal conductivity, x, 
were selected by taking the turbulent Schmidt, Sc, and Prandtl, Pr, numbers 
to be equal to unity (refs. 14 and 15) 

yC 

Sc = 1 Pr = 1 (35) 

pO^B 

One set of calculations were performed in which the effective transport prop- 
erties (y and x) were reduced by a factor of 100. The results obtained with 
the different transport properties (but with the same geometric and operating 
conditions) showed the flow patterns for the two problems to be almost 
identical. 

The axial and radial components of the gas velocities inside the cylinder 
were calculated as functions of time and represented graphically. Only the 
graphical results were presented here. The resultant velocities are represen- 
ted by arrows, drawn to scale when the magnitude of the velocity vectors are 
less than or equal to 0.5 meter per second. When the magnitude of the velocity 
vector is greater than 0.5 meter per second, it is set equal to 0.5 meter per 
second. 


Single Component Gas - Intake Stroke 

The results are shown in figures 2 to 6. Throughout the intake stroke, a 
jet issued from the valve opening. For the conditions employed in this study, 
the jet always impinged on the piston surface (figs. 2 to 6). The existence 
of such a jet has been observed previously in both experimental (refs. 16 to 
22) and numerical (refs. 23 to 38) studies of similar piston-cylinder problems. 

During the first part of the intake stroke (crank angle i) less than 
~10 degrees measured from the TDC position), only the jet existed inside the 
cylinder and there were no recirculating flows. The reason for this is that 
the pressure across the valve opening was taken to be constant. This allowed 
the jet to spread out radially from the valve opening without causing the jet 
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to separate. In an actual piston-cylinder, where the static pressure across 
the valve opening would not be constant, recirculating flows may be established 
earlier than is accounted for in this study. 

At some crank angles greater than ~10 degrees (but always less than 90 deg- 
rees), flow separation produced two recirculating flows having the shape of 
toroidal vortices. One of these vortices (referred to as the cylinder-head 
vortex) was located between the jet and the cylinder wall. The other (referred 
to as valve vortex) was located between the jet and the center line. The for- 
mation of such vortices during the intake stroke has been observed previously 
both in experimental (refs. 17 to 22) and numerical (refs. 23 to 39) studies. 

The formation and subsequent behavior of these vortices depended largely 
on the speed and angle of the jet issuing from the valve opening. The jet 
speed is a function of the pressure difference across the valve opening (i.e., 
the difference between the constant pressure of the intake charge P-j and the 
instantaneous pressure of the gas inside the cylinder). On one hand, the pres- 
sure inside the cylinder increases with time due to the influx of the gas. On 
the other, the pressure inside the cylinder decreases with time due to the mo- 
tion of the piston away from the cylinder head. 

When the intake charge pressure was the same as the residual gas pressure 
(P-j/Pc = 1), the jet speed continuously increased as the piston moved away 
from the cylinder head until the crank angle jJ reached 90 degrees. Thus, 
once the cylinder-head and valve vortices were formed, the jet continuously 
strengthened these vortices, and the vortices steadily grew in size through- 
out the intake stroke. The growth of the vortices during the intake stroke 
was observed previously in experimental (refs. 17 to 21) and numerical (refs. 

26 to 38) studies. 

Previously, Ekchian and Hoult (ref. 17) reported the speed with which the 
cylinder-head and valve vortices moved. These investigators performed exper- 
iments with a valve-seat angle of 45 degrees, a piston-to-valve opening area 
ratio of 23.5, bore-to-stroke ratios of 0.75 to 1.35, and crankshaft angular 
speeds of 750 to 2500 rpm. Ekchian and Hoult noted that: (a) the valve vor- 

tex grew in size and its center of rotation moved at about one-half of the 
piston speed during the intake stroke; (b) the cylinder-head vortex was ap- 
proximately constant in size and was almost stationary during the intake 
stroke; and (c) the cylinder-head vortex was trapped in the volume formed by 
the jet, the cylinder head, and the cylinder wall. 

Ekchian and Hoult' s results can be compared with the numerical calcula- 
tions performed here for valve seat angle a equal to 45 degrees, piston-to- 
valve-opening area ratio of 4, bore-to-stroke ratio of 1.67, and crank-shaft 
angular speed of 400 rpm. These numerical results show the following phenom- 
ena: (a) the valve vortex grew in size and its center of rotation moved at 
about one-third to one-half of the piston speed as the piston moved away 
from the cylinder head (fig. 2). This is in fair agreement with the experi- 
mental observations of Ekchian and Hoult (ref. 17), and (b) the cylinder-head 
vortex grew in size and its center of rotation moved at about one-fifth of 
the piston speed (fig. 2). 

The growth of the cylinder-head vortex and its corresponding motion re- 
sulted from the fact that the jet bent towards the piston. This bending 
increased the volume available for the vortex to occupy. In Ekchian and 
Hoult' s experiment, the jet remained straight; hence the size of the 
cylinder-head vortex remained constant and its center of rotation did not 
move. The reason that the jet did not bend in their experiment is probably 
due to the fact that the jet speed in their experiment was higher than the 
jet speed simulated in this study. 
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When the intake charge pressure is higher than the residual gas pressure 
(Pi/Pq > 1), then, depending of the piston speed, the speed of the jet 
may first increase, then decrease, and then increase again. During the first 
jet-speed increase, both the cylinder-head and valve vortices formed. De- 
pending on the jet speed and the strength of the cylinder-head and valve 
vortices formed during the first jet-speed increase, these vortices may 
diminish in strength or even completely die out during the subsequent de- 
crease in jet speed. When the vortices die out, a new set cf cylinder-head 
and valve vortices form when the jet speed increases again due to the piston 
motion (fig. 3). It is interesting to note that the flow patterns inside 
the cylinder formed with the intake-charge-to-residual-gas pressure ratios 
of unity and 1.036 become similar at crank angles beyond ~40 degrees (figs. 

3 and 4) . 

The cylinder-head and valve vortices formed at different times. At 
higher piston speeds (and corresponding higher jet speeds), both the 
cylinder-head and valve vortices formed earlier. When the valve-seat angle 
a equals zero degrees, the cylinder-head vortex always formed before the 
valve vortex (figs. 4 to 6). When the valve-seat angle a equals 45 deg- 
rees, the valve vortex formed before the cylinder-head vortex (fig. 2). 

The relative strengths of the cylinder-head and valve vortices depended 
upon the relative times at which these vortices form. When the valve vortex 
formed at crank angles near 90 degrees (i.e., at the point at which the jet 
speed began to decrease due to decreasing piston speed), the cylinder-head 
vortex was always stronger than the valve vortex (figs. 3 and 4). When the 
valve vortex formed early (0 = 30 to 60 degrees), the valve vortex became 
stronger than the cylinder-head vortex by the end of the intake stroke (figs. 
2, 5, and 6) . 

At a low piston speed (n = 400 rpm and rc = 0.03 meter) a third 
recirculating flow, in the shape of a toroidal vortex, formed between the 
valve vortex and the valve surface (fig. 4). Such a vortex was observed 
previously at low piston speeds both in experimental (refs. 19 and 21) and 
in numerical (refs. 37 and 38) studies. 

At 400 rpm, bore-to-stroke ratio of 0.78, and valve-seat angle of zero 
degrees, the cylinder-head vortex broke up into two smaller toroidal vortices 
enclosed by a larger toroidal vortex (fig. 6). Ekchian and Moult (ref. 17) 
observed the break-up of the valve vortex into a random flow. In the present 
study, the vortex did not degenerate into a random flow because instability 
and three-dimensional flow (needed for the complete vortex break-up) were not 
included in the numerical solution. 


Single Component Gas - Compression Stroke 

During the compression stroke, the cylinder-head and valve vortices were 
no longer separated by the jet, allowing the two vortices to interact (figs. 

7 to 10). Because of the rotational motion of the two vortices, they forced 
each other towards the piston. Since the two vortices were of unequal 
strength, the weaker one was pushed closer to the piston surface by the 
stronger one. 

As the piston speed increased during the compression stroke (crank angle 
i between ~190 and 240 degrees), the two vortices began to coalesce. 

Coalesce as used here means combine by entrainment. The time of this coales- 
cing depended upon the strength of the vortex which was nearer to the piston 
surface. 
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At crank angles 0 between ~230 and 240 degrees, the two vortices 
coalesced into a single toroidal vortex. The coalescing of the cylinder-head 
and valve vortices has not been described by previous investigators. How- 
ever, the disappearance of one vortex was found by Ashurst (ref. 26) and by 
Oiwaker, et al. (ref. 27) in numerical simulations of similar piston-cylinder 
problems. 

Towards the end of the compression stroke (crank angle ii > 320 deg- 
rees) a new recirculating flow (corner vortex) formed in the corner of the 
cylinder head and cylinder wall. The numerical simulations of Chong, et al. 
(ref. 40) and Gosman, et al. (ref. 37) also showed the formation of the cor- 
ner vortex. Both this study and Chong, et al. (ref. 40) indicate that the 
corner vortex formed earlier at lower piston speeds. 

In an actual cylinder, a new vortex forms near the cylinder wall due to 
the piston scraping off the boundary layer next to the cylinder wall during 
the compression stroke (refs. 1 and 41 to 46). This vortex was not observed 
in the present study. One reason for this is that a large number of grid 
points would be required for a high enough resolution to observe this 
phenomenon. 


Two-Component Gas Mixture 

Numerical solutions were also obtained when the intake charge is a two- 
component gas mixture of components A and B. The mass fractions of compo- 
nents A and B at the valve opening varied with time according to the expres- 
sion given in table 2. The mass fraction are represented graphically by 
symbols selected to give reasonable visual contrast for the different mass 
fractions inside the cylinder. The darker and lighter regions represent 
regions of lower and higher X/\, respectively. 

Whether the intake charge was a single-component gas or a two-component 
gas mixture, the flow pattern inside the cylinder depended largely on the 
speed of the jet issuing from the valve opening during the intake stroke. 

In single-component gas problems, the density of the intake charge was con- 
stant and the jet speed depended only on the pressure difference across the 
valve opening. In two-component gas mixture problems, the intake charge 
density varied with time so that the jet speed depended not only on the 
pressure difference across the valve opening but also on the intake charge 
density. 

As in single-component gas problems, a cylinder-head and a valve vortex 
formed. Since the jet speed varied with time during the intake stroke (as 
it did in single-component gas problems), the vortices formed during the 
early part of the intake stroke may weaken and then strengthen again, or may 
die out completely and then form again (fig. 11). During the latter part of 
the intake stroke and during all of the compression stroke, the essential 
features of the flow were similar for single-component and two-component gas 
mixtures (fig. 12). 

At the beginning of the intake stroke, the gas inside the cylinder (the 
residual gas) was a single-component gas with the same properties as compo- 
nent A of the intake charge. As the intake charge (composed of components A 
and B) entered into the cylinder, mixing of the gases took place. As time 
progressed, the composition of the gas mixture inside the cylinder tended to 
become more uniform due to mass diffusion and convective flow. The degree of 
mixing depended on the spatial gradient of the mass fraction of component A, 
the speed with which the gas was convected to the various regions inside the 
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cylinder, and the time available for mass diffusion and convective flow. 
These parameters in turn depended on the geometry, the crankshaft angular 
speed, the intake-charge-to-residual-gas pressure ratio, and the mass frac- 
tion of component A at the valve opening (i.e., the nonhomogeneity of the 
intake charge). For the conditions employed in this study, the gas mixture 
inside the cylinder always became "homogeneous" before the end of the com- 
pression stroke (figs. 13 and 14). By "homogeneous," it is meant that the 
mass fraction of component A inside the cylinder was uniform to within plus 
or minus of one percent. 


CONCLUDING REMARKS 

The results generated in this study illustrate some of the major features 
of the flow inside an axisymmetric piston-cylinder configuration. The re- 
sults were obtained by equations formulated for axisymmetric, laminar flows. 
From the results obtained thus far, it appears that it would be worthwhile 
to extend the implicit-factored method to the study of turbulent flows. 
Extension of the method to three-dimensional internal flows also appears 
worth exploring. Application of this method to three-dimensional piston- 
cylinder problems could shed light on the three-dimensional hydrodynamic 
instabilities and on the vortex break-up. 
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TABLE I. - GOVERNING EQUATIONS* 



*Cp-j = constant pressure specific heat of component i, D/\g = binary diffusion coeffi- 
cient for components A and B, e = energy per unit volume, h-j = specific enthalpy of 
component i, h^ = standard enthalpy of formation per unit mass for component i at 
temperature Ti, = molecular weight of component i, P = static pressure, qj = heat 
flux in the j-directi^on, R = universal gas constant, Vj = j-component of the velocity, 
V^-j = j-component of ''the diffusion velocity for component i, = mass fraction of 
component i, x = thermal conductivity, y = viscosity, p = density, t = shear stress. 


12 






TABLE I. - Concluded. 



^66 = 2u/-|y(7 • v) 


^zz’ ir-f 


^rz " 1 3z 3r 


1 3 ‘*''z 

V . V = - — rV + — ^ 
r 3r r 3z 


Xb = 1 - Xa 


^ ^A*^A * ^B*^B ^ ^A^*^A “ ^B^ ^ *^B 


'’a ■ ''A <^P^ - '’A ^ <:pA(T - Ti) 


Equation 

number 


he = h 




CpB dT = hg + Cpg(T - T^) 


* - X ^ ^ ('’a ■* ^’b) ^^PA ■ ~ ir ^A 

~ ^ 3? ‘’Vaz^'^A “ “^B^ " “ ^ It “ ^ (“^A ~ “^b) ^^PA " It h 











TABLE 2. - BOUNDARY AND INITIAL CONDITIONS 


Equation 


V = O’ 


V, ^ = 0 


T = T. 


T = T. 


T = T, 


w 


T = T, 


^ X„ = 0 


3Z ''A 


ar "A 


X„ = 0 


P= P^, T= T. 


Xa = Xa <0 


(l - j s1t> ( 


20 


Xfi = 1 - Xa 


^ - \J tan a 
r z 


le. + i i_ r(pV ) + — (pV ) = 0 
at r ar az z' 


aV. 


ie. w "’z ae Q 

ar " ar " r " ar " ar 


p = P(,/(RT./M^), X^ = 1, = 


^ ■ fc [<¥pa - ^ W - =paTi 



Equation 

number 


Piston surface 

20 


Cylinder head, cylinder 
wall, and valve 

21 


Valve 

22 


Cylinder head 

23 


Cylinder wall 

24 


Piston surface 

25 


Cylinder head, valve, 
and piston surface 

26 


Cylinder wall 

27 


Valve opening 

28 


Valve opening 

29 


Valve opening 

30 

j 

Valve opening 

31 


Valve opening 

32 


Center line 

33 

0 

)] 

(Everywhere inside the 
J piston-cylinder con- 
1 figuration at time 
t = 0 

34 
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TABLE 3. - SUMMARY OF PARAMETERS USED IN THE NUMERICAL CALCULATIONS* 


Type of 
problem 

p c 

n(rpm) 


a(deg) 

o)(in eq. (29)) 

I^(in eq. (29)) 

One- 

1.67 

400 

1 

45 

0 

1 

component 

1.67 

400 

1.036 

0 

0 

1 

gas 

1.67 

400 

1 

0 

0 

1 


1.67 

1000 

1 

0 

0 

1 


.78 

400 

1 

0 

0 

1 

Two- 

1.67 

1000 

1 

0 

1 

0.938 

component 

1.67 

400 

1 

0 

0 

.938 

gas 

1.67 

400 

1 

0 

1 

.938 

mixture 

.78 

400 

1 

0 

1 

.938 


1.67 

400 

1.036 

0 

1 

.938 


3 

*The following parameters were the same for all cases studied: u = 3.2xl0~ 

kg/m-s, X = 3.224 kg-m/s^-K, h° = 42,252 m^/s^, h° = 1.7x10^ m^/s^, 

Cp^ = 1,006.6 m^/s^-K, Cpg = 1.653.6 m^/s^-K, M^ = 28.96 kg/kg-mole. 

Mg = 114.2 kg/kg-mole, R = 8314.3 kg-m^/kg-mole-K-s, D^g = 3.2xl0”^/p 

m^/s^, r_ a 0.05 m, r , = 0.01875 m, r^ = 0.03438 m, 6(t = 0) = 0.01 m, 

p V n 

0-2 m, T. = T^ = T^ a T^ a Tp = 340 K, and a 101,325 kg/m-s^. 
Results are shown only for the first six cases listed in this table. 
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Figure 5. - Flew patterns during the Intake stroke as a function of 
crank angle from TOC for single-component gas, ■ L67, 
n ■ 1000 rpm, Pj/Pj - 1, and a • 0° (Table 3). 
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Figure 6. - Flow psttems during the intake stroke as a function of crank angle from TDC 
for single-component gas, fj/rj • 0.78, n • 400 rpm, P|fPc • 1, anda ■ oP (Table 3). 


Figure 7. - Flow patterns during the compression stroke as a function 
of crank angle from TOC for single-component gas, rJrj- L67, 
n • 400 rpm, P|fPp ■ L036, and a • 0° (Table 3). 
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Figure 10. - Flow patterns during ttie compression stroke as a function of crank angle 
from TDC for single-component gas, rJrj. ■ 0.78, fJ ■ 400 rpm, Pj/P^ • 1, ando ■ 0® 
(Table 3). 



Figure IL - Flow patterns during the Intake stroke as a function of 
crafik angle for two-component gas mixture, u • 1, Urj • L67, 
0 ■ 1000 rpm, P|/Pj. • LO, ando ■ oP (Tables). ^ 
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Figure 13. - Distribution of air (component A) to fuel (component 81 
ratios during ttie Intake stroke as a function of crank angle from 
TOC for two-component gas mixture, u • 1. rJrj ■ L67, fl ■ 1000 
rpm, Pj(Pp ■ LO, ando ■ 0® (Table 3). Darker regions Indicate 
lower alr-to-fuel ratios. 
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Figure 14. - Distribution of air (component A) to fuel (component B) 
ratios during the compression stroke as a function of crank angle 
fromTDC for two-component gas mlxture.w ■ 1, rJtr ‘ 1.67, 
n ■ 1000 rpm, P|/Pg ■ LO, anda • oP ITable3.) Darker regions 
Indicate lower alr-to-fuel ratios. 
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